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Abstract 

We consider the scenario in which the mass matrices of the charged fermions 
in the SO(10) Grand Unified Theory are generated exclusively by renormalizable 
Yukawa couplings to one 10 © 120 representation of scalars. We analyze, partly 
analytically and partly numerically, this scenario in the three-generations case. We 
demonstrate that it leads to unification of the b and r masses at the GUT scale. 
Testing this scenario against the mass values at the GUT scale, obtained from the 
renormalization-group evolution in the minimal SUSY extension of the Standard 
Model, we find that it is not viable: either the down-quark mass or the top-quark 
mass must be unrealistically low. If we include the CKM mixing angles in the test, 
then, in order that the mixing angles are well reproduced, either the top-quark mass 
or the strange-quark mass together with the down-quark mass must be very low. We 
conclude that, assuming a SUSY 50(10) scenario, charged-fermion mass generation 
based exclusively on one 10 © 120 representation of scalars is in contradiction with 
experiment. 



*E-mail: balio@cftp.ist.utl.pt 
^E-mail: helmut.kuehboeck@gmx.at 
■^E-mail: waiter. gr imus@ uni vie . ac . at 



1 



1 Introduction 



All the fermionic multiplets of one Standard-Model generation, plus one right-handed neu- 
trino singlet, fit exactly into the 16-dimensional irreducible representation of the Grand- 
Unification group 5*0(10). This is the unique and distinguishing feature of the unified 
gauge theories (GUTs) based on this group [I]. As a bonus, the presence of three (for 
three generations, as we shall assume in this paper) right-handed neutrino singlets allows 
one to incorporate into this GUT the seesaw mechanism of type I [2] . 

However, when it comes to the scalar sector and to fermion mass generation the 
uniqueness of the 50(10) GUT is lost and numerous ramifications exist. One possible 
strategy to limit the freedom in the scalar sector is to confine oneself to renormalizable 
terms — for a review see, for instance, In that case, the scalar representations coupling 
to the fermions are determined by the relation jH U] 

16® 16= (10©126) s ffil20 AS , (1) 

where the subscripts "S" and "AS" denote, respectively, the symmetric and antisymmetric 
parts of the tensor product. Thus, scalars with renormalizable Yukawa couplings to the 
fermions must transform under 50(10) either as 10, 126, or 120 (the 10 and 120 are 
real representations; the 126 is complex). A minimal supersymmetric (SUSY) scenario — 
which has built-in the gauge-coupling unification of the minimal SUSY extension of the 
Standard Model (MSSM) — making use of one 10 and one 126 for the Yukawa couplings [6J 
has recently received a lot of attention. This attention was triggered by the observation [7j 
that maximal atmospheric neutrino mixing may in this theory be related to b-r unification 
via the type II seesaw mechanism [SJ. Detailed and elaborate studies of this minimal 
theory have been performed for its Yukawa couplings (11110111111121 and scalar potential O 
IT%j . This "minimal SUSY 5O(10) GUT" works very well, since its Yukawa couplings are 
able to fit all fermion masses and mixings, allowing in particular for small quark mixings 
simultaneously with large leptonic mixings. However, in this context a minimal Higgs 
scalar sector is too constrained jT5] and does not allow to produce large enough neutrino 
masses [16]. 

As a way out, the 120 scalar representation — which had been somewhat arbitrarily 
left out — may be used for a rescue ^3 QB] • in OH EE] that representation was only taken 
as a perturbation of the minimal scenario, to cure minor deficiencies in the fermionic 
sector. However, in ^1 it was pointed out that the antisymmetric coupling matrix of 
the 120 could be responsible for the different features of quark and lepton mixing, since 
that matrix has different weights in all four Dirac-type mass matrices — i.e. in the Dirac 
mass matrices for the up-type quarks, down-type quarks, charged leptons, and neutrinos. 
Thinking along this line, the roles of the 120 and 126 could be interchanged in the 
charged-fermion sector: the brunt could be borne by one 10 © 120, and the Yukawa 
couplings of the 126 would be just a perturbation. This thought is realized in the model 
of ^7j, where the scalar 126 is still a protagonist in the neutrino sector, through the 
type I seesaw mechanism. 

In this paper we investigate the extreme form of this scenario of ^Zj, namely we assume 
that the 126 plays no role whatsoever in the Yukawa couplings to the charged fermions, 1 

1 This idea was previously put forward in |2(J) . 
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and may be important only in the neutrino sector, where it would responsible for large 
Majorana neutrino masses. Thus, we base our investigation on the following assumptions: 

i) The charged-fermion mass matrices result solely from the Yukawa couplings of one 
10 and one 120 scalar multiplets. 

ii) The mechanism for the generation of the light-neutrino mass matrix is the type I 
seesaw mechanism, possibly with some admixture of type II. 

Due to the first assumption, the Yukawa-coupling matrix of the 126 can be used freely for 
the neutrino mass matrix and, therefore, one can accommodate any neutrino masses and 
lepton mixing that one wants, through either the type I or type II seesaw mechanisms. 
The tight connection between the charged-fermion and neutrino sectors is lost, and the 
predictive power of the model for the neutrino sector too. The subject of this paper is 
then only the discussion of the charged-fermion masses and of quark (CKM) mixing under 
the assumption i), and the working out of where this assumption is successful and where 
it might fail. 

The charged-fermion sector and the Yukawa couplings in [IS] coincide with ours. Still, 
our results do not, in general, apply to that model. The reason is that its authors assume 
split supersymmetry, where the renormalizat ion-group evolution of the fermion masses 
differs from the one of the MSSM. Indeed, in order to test any specific scenario one must 
use the charged-fermion masses and the quark mixings at the GUT scale. Having in mind 
a SUSY 5*0(10) GUT and the MSSM, we use in this paper the values computed in [2T] 
with the renormalizat ion-group evolution of the MSSM. 

This paper is organized as follows. In Section El we discuss the mass matrices and 
count the number of parameters. Basis-invariant quantities are introduced in Section 01 
The derivation of some inequalities, and b-r unification, are discussed in Section 0] In 
Section El we show that a partly analytical treatment of our scenario is possible when the 
Yukawa-coupling matrices are assumed to be real. Section El explains our procedure for 
the numerical fit of the mass matrices to the fermion masses and to the CKM mixing 
angles at the GUT scale. We present our results in Sectional which is followed by a brief 
summary in Section |H1 

2 The charged-fermion mass matrices 

The mass Lagrangian that we are concerned with is 



The symmetric and antisymmetric Yukawa couplings of one 10 and one 120 scalar repre- 
sentations, respectively [3], generate the mass matrices, which at the GUT scale may be 
parametrized as 



C M = -d L M d d R - I l MiIr ~ u L M u u R + H.c. 



(2) 



M d 



S + e^A, 
S + re 19 A, 
P S + qe*A, 



(3) 
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S being symmetric while A is antisymmetric. The parameters p, q, and r are real and 
positive. The matrix S is proportional to the Yukawa-coupling matrix of the 10, while 
A is proportional to the Yukawa couplings of the 120. The factors e 1 ^, re t6 , p, and qe 1 ^ 
depend on some ratios of vacuum expectation values. 
We may perform changes of weak basis 

S US IF, 

A _ UAU T , { ' 

where U is unitary. In this way we may reach convenient weak bases. We may for instance 
use U to diagonalize S: 



(5) 



with real and non-negative a, b, and c. Alternatively, we may use U to force A to have 
only two non-zero matrix elements, and moreover two matrix elements of S to vanish: 
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(6) 



In the weak basis (jfi|L again, we may choose a, b, and c to be real and non-negative. 

As for the number of degrees of freedom in the mass matrices 0, we consider two 
cases: 

• In the complex- Yukawa-couplings case, the matrices S and A are a priori complex 
and contain nine independent matrix elements (six in S and three in A), hence nine 
moduli and nine phases. One of three phases ip, 8, and £ may be absorbed in the 
definition of A. Through a weak-basis transformation we may eliminate the three 
moduli and six phases which parametrize U. In that case the model has, therefore, 
nine real parameters and five phases. 

• In the real- Yukawa-couplings case, in which CP violation is considered to be sponta- 
neous, the matrices S and A are a priori real and contain nine independent moduli. 
If we want to preserve the reality of S and A, the matrix U of the weak-basis 
transformation Q must be chosen real (orthogonal), 2 hence it contains three real 
parameters. One ends up with nine real parameters as before, but only three phases. 

With these 14 (in the complex case) or 12 (in the real case) parameters we must try and 
fit 13 observables: nine charged-fermion masses and four parameters of the CKM matrix. 
Even if there is, in the complex case, an excessive number of parameters, the fitting may 
prove impossible, due to the fact that a large number of those parameters are phases. 

2 If S and U are assumed to be real, then one may obtain a weak basis of the form JSJ, but a, b, and c 
must be allowed to be negative. It is only when we allow U to include some i factors that we may obtain 
non-negative a, b, and c; but then x, y, and z will not necessarily be real. 
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3 Fermion masses and invariants 

We first confine ourselves to the masses. For brevity of notation we introduce 

Pd 



m 2 d + m 2 s + ml, 



m 2 d m 2 s + m 2 ml + m 2 m\, 



ml + m 2 + m 2 , 
m 2 + ml + ml, 



m l m fi + rfjjn 2 - + m 2 T m 2 e , 
m 2 u m 2 c + m 2 c ml + m 2 t m 2 u , 



ltd 



m 2 d m 2 s m 2 b , 

111 
mjnjn T , 

111 
m u m c m t . 



(7) 



We define the matrices H a = M a M\ (a = d, £, u), which have eigenvalue equations 

det (m 2 t — H^j = m 6 — <j a m 4 + p a m 2 — ir a = 0. 
With the mass matrices we obtain the relations 

a d = s 2 + 2a 2 , 
oi = s 2 + 2r 2 a 2 , 
<y u = P 2 s 2 + 2q 2 a 2 , 



where 



also, 



where 



finally, 



where 



Pu 



s 2 = tr {SS*), 
a 2 = -~tv(AA*); 



s 4 + a\ + 2z 4 + 2 Re (e 2i % 4 ) , 

s 4 + r A a 2 2 + 2r 2 z A + 2r 2 Re (e 2ie z 4 ) , 

p 4 s 4 + q A a 2 2 + 2p 2 q 2 z A + 2p 2 q 2 Re (e 2 % 4 ) 



1 



s 4 

2 4 



s 2 2 - tr (SS*SS*) 



7T,, 



s 2 a 2 + tr (SS*AA*) 
~ti (AS* AS*); 



S3 + e 2 ^z 3 
s 3 + r 2 e 2W z 3 



p 3 s 3 + pq 2 e 2i ^z 3 



S3 

z 3 



det S, 



tr (SA 2 ) -ItrStr (A 2 ) 



(9) 
(10) 

(ID 



(12) 
(13) 

(14) 
(15) 
(16) 

(17) 
(18) 
(19) 

(20) 
(21) 
(22) 

(23) 
(24) 
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4 b—r unification 



In pH] the mass matrices for the charged-fermion sector are the same as in this paper, but 
the discussion is confined to the two-generations case. In that paper, approximate b-r 
unification is traced back to some inequalities derived from the specific structure of the 
mass matrices. Here we show that analogous inequalities hold in the three-generations 
case. 

It is convenient to use the weak basis of Eq. (jSJ). We remind that, in that weak basis, 
a, b, and c are real and non-negative, while x, y, and z are in general complex. One has 



S-2 

a 2 
s 4 
za 

Z4 
S3 
Z3 



2 i 1.2 i 2 

a +o + c , 

I 1 2 . I i2 . I 1 2 

\x\ + \y\ + \z\ , 

2 iS2 , 1.2 2 , 2 2 

a o + be + c a , 

2 112,72112, 2 I I 

a \x\ +o \y\ + c \z\ 
bcx 2 + cay 2 + abz 2 , 
abc, 

2 i i, 2 , 2 

ax + by + cz . 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 



Note that z^ and z^ are in general complex, while the other parameters are real. 
From Eqs. ^B>, flU), and fl2EJ) we derive 



Pu > p 4 S4 + q^a\ + 2p 2 q 2 Z4 — 2p 2 q 2 (be \x\ 2 + ca \y\ 2 + ab \z 

= p 4 s 4 + q^a\ + 2p 2 q 2 (a 2 — bej \x\ 2 + (b 2 — caj \y\ 2 + (c 2 — ab) \z\ 
Without loss of generality we assume that 

b> a, c> a. 

Since a + b + c is non-negative, the inequalities (jHlj) are equivalent to 

b 2 — ca > a 2 — be, c 2 — ab > a 2 — be. 
Applying the inequalities (|33|) to the inequality (}3*3*|) and remembering Eq. 

Pu > p i SA + q^a\ + 2p 2 q 2 (a 2 — bej a 2 . 
We next rewrite Eqs. (fTTj) and (|23j) as 

° 2 = 2g2 au ~ p2 (° 2 + &2 + ° 2 ) ■ 
We plug this equation into inequality (|3l)j) and find after some algebra that 



(32) 
(33) 

(34) 



where 



Pu>~ 



p 2 a 2 



p 2 (b + cf + F, 



o u +P 2 (b + c) 



2^2 



3p z a 



(35) 
we obtain 
(36) 

(37) 

(38) 
(39) 
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The inequalities give b + c > 2a, hence 

F >^(a u + Xf) (40) 

> t**L (41) 

> 0. (42) 

From inequalities (jHEJ) and 

Pn> J [a u -p 2 (6 + c) 2 ] 2 . (43) 
This inequality may equivalently be written 

<p 2 {b + cf <a u + 1Jp- u . (44) 
It is obvious that, in an exactly analogous fashion, one may derive 

(45) 



(Td-Zy/pd < (b + c) < a d + 2^fp~ d , 
a e - 2y/pi < (b + c) 2 < a £ + 



Inequalities (|45|) should be compared with those of ^H] ■ One reaches the same conclusion 
as in [THj: the intervals [a d — 2^fp~ d , a d + 2^fp~ d \ and [at — 2J~pl, at + 2^/pt] must overlap. 
This overlap — at the GUT scale — implies that, at that scale, mj ~ m T . Notice that this 
conclusion was reached without making use of the quantities n a . 

Comparing inequalities (J44j) and (J45|) . one also finds that the parameter p is approxi- 
mately given by 

p ~ — (46) 
m b 

at the GUT scale. 

Inequality (|38p also delivers F < p u . Taking into account inequality (j41|) . one has 

P 2 a 2 < ^. (47) 
With Eq. pfijl in mind, this gives, approximately, 

a < (48) 
m t 

Numerically, using the values of the quark masses in the MSSM at the GUT scale, as 
given in [21], one obtains for instance a < 3.8 Me V for tan/5 = 10. 
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5 Analytical treatment of the real case 

In this section we analyze the case of real Yukawa-coupling matrices, i.e. the case of real 
S and A. In this case it is convenient to define 



%2 



Then, 



«2 



tr S, 

1 r 2 



tr 



tr S 



(49) 
(50) 



( S 2 ) 



x x — 2x 2 



S4 — — 2,XiS^. 

With 5 and A real, 24 is real, and moreover it is not independent from z 4 , rather 

Z4 — Z 4 = X2CL2 ~ x l z 3- 

This allows one to write Equations jHj) and (|T5|) as 



(51) 
(52) 



(53) 



ft 



ft 



x 2 — 2X1S3 + a 2 , + 2Z4 + 2 cos (2ijj)z 4 

^2 — 2xis 3 + a 2 — 2x 2 a 2 + 2x x z 3 + 2 [1 + cos (2-0)] z 4 

(x 2 - a 2 ) 2 + 2xi (z 3 - s 3 ) + 2 [1 + cos (2ip)] z 4 , 

x\ - 2x1 s 3 + r A a\ + 2r 2 ^ 4 + 2r 2 cos (26»)z 4 

x\ - 2xis 3 + r 4 a2 - 2r 2 x 2 a 2 + 2r 2 XiZ 3 + 2r 2 [1 + cos (26)} z 4 

x 2 - r 2 a 2 ) 2 + 2xi (r 2 z 3 - s 3 ) + 2r 2 [1 + cos (26)] z 4 . 



Now, plugging Equation (J5T|) into Equations (jHJ) and (jTUjl . one obtains 

cr d = x 2 - 2 (x 2 - a 2 ) , 
= x 2 — 2 (x 2 — r 2 a 2 

Hence, Equations ()54|) and f)55|) may be rewritten as 



Pd 



1 ( 2 
7 (fx ~ vt) 



2xi (23 - s 3 ) + 2 [1 + cos (2^)] z 4 , 
2xi (r 2 -2 3 - s 3 ) + 2r 2 [1 + cos (26)] z A . 



(54) 



(55) 



(56) 
(57) 



(58) 
(59) 



In the trivial case cos (2ip) = cos (26) = —1, the mass matrices and Mg are Hermitian 
and their eigenvalues directly yield the fermion masses. Discarding that rather trivial case 
from consideration, we find that Equations (}58j) and (}59|) lead to 



= r 2 [1 + cos (26)} 
- [1 + cos (2V>)] 



Pd 



x i - a d) + 2xi (s 3 

2 



^3 



- - (x\ - a e ) + 2a;i (s 3 - r 2 z^ 



(60) 



On the other hand, since S3 and z 3 are real when the matrices S and A are real, 
Equations (|20|) and (|2*TJl read in that case 



S l + Z t + 2 COS (2l^)s 3 Z 3 = TTd, 

s 3 + r 2 3 + 2r cos (2^)53^3 = 7T£. 



(61) 



Defining 

/1 = 1-r 4 , (62) 

/ 2 = cos (2^) -r 2 cos (20), (63) 

f 3 = r 4 cos (2^>) - r 2 cos (20) , (64) 

/ 4 = 7Td-VT£, (65) 

/5 = rVrf - 7T£, (66) 

/ 6 = r 2 cos (20)^ -cos (2^ K (67) 
the system of equations (}6T|) has solutions given by 



S 3 



■/l/5 _ 2/ 3 / 6 ± 2/ 3 y /g — / 4 /r 

A 2 + 4/2/3 



2 /1/4 — 2/2/6 T 2/2V /e — AA 

23 = JfTIM • (68) 



S3Z3 



A A + A A =•= Ay A ~ A A 

A 2 + 4AA 



We use as input the three charged-lepton masses, the three down-type-quark masses, 
and also r, cos (29), and cos (2?/>). Equations allow us to compute S3 and z 3 from that 
input. Inserting those values of S3 and z% in Equation (JBTij) . we obtain a quartic equation 
for xi, which may be analytically solved. The quantities S2 and a 2 are then computed as 

x\ r 2 a d - a t 

32 = Y + W^-y (69) 

fl2 = 27J^)- (70) 

Finally, z 4 is computed from either Equation or Equation (J59*j) . and z 4 is obtained 
from Equation (|53|). All the invariants pertaining to the matrices S 1 and A are thus 
analytically computed from the input. 

One must, yet, take into account the fact that those invariants must satisfy several 
inequalities. In the weak basis (jSJ, 

x\ = a + b + c, (71) 
x 2 = ab + be + ca, (72) 
S3 = abc. (73) 
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The numbers a, b, and c are real, and they may be negative, see footnote 121 The quantity 



A = x\x\ + 18xix 2 s 3 - - 4a; i s 3 - 27S3 (74) 
= \{a-b)(b-c){c-a)f (75) 

must therefore be non-negative. Further non-negative quantities may be conveniently 
derived by using the weak basis (jOJ and deriving, in that basis, the values of a, b, c, d 2 , 
and f 2 from the invariants. From the condition that f 2 must be non-negative one obtains 

£ = a 2 z^ — z 2 (76) 
> 0. (77) 

From the condition that d 2 must be non-negative one obtains 

v]/ = -z\ + z\ (2xiz 3 - x 2 a 2 ) + z A a 2 z 3 (3s 3 + Xix 2 ) - Xi-s 3 a 2 2 - (x\ + x 2 ) z\ 

+zf (xix 2 - s 3 ) - a 2 z\ (x\ + X1S3] + 2x 2 s 3 alz 3 - sla\ (78) 
> 0. (79) 

The conditions that A, E, and ^ be non- negative constitute a severe constraint on the 
inputted values of the charged-fermion masses and of r, 6, ip. 

After having computed the invariants, one may further input the three up-type-quark 
masses and therefrom derive the values of p 2 , q 2 , and cos (2£). In practice, this involves 
solving a cubic equation, and thereafter imposing the constraints p 2 > 0, q 2 > 0, and 
I cos (2£)| < 1. This obviously translates into constraints on the inputted up-type-quark 
masses. 

In the way delineated in this section, one may analytically solve the case of real S and 
A matrices, by inputting the charged-fermion masses and therefrom deriving S and A, 
without having to have recourse to fits. In practice, however, doing things the other way 
round — trying to fit the charged-fermion masses numerically from some inputted values 
of S, A, and the other parameters — proves more effective. We turn to that procedure in 
the next section. 



6 The fitting procedure 

In order to check whether the mass matrices allow to reproduce the masses and CKM 
mixing angles at the GUT scale, we use a x 2 analysis, as was previously applied for 
instance in (2HIE1]- As for the masses, the x 2 -function is given by 

Xmasses = xl + xl + xl, (80) 

where 

2 x - ( mdx) - fhi 



xi= £ • (81) 

i=d,s,b \ / 

and analogously for xj u - The masses at the GUT scale are fhi ± Srrii, whereas the 
rrii(x) are the masses calculated from Eqs. (j2J) as functions of the parameter set x = 



10 



{S, A,p, q, r, if), 9, £} (see Section |21 for the distinction between the "real" and the "com- 
plex" cases). The total x 2 -function is the sum 



with 



Xtotal — Xmasses + XcKM' (§2) 

2 x ^ ( sin 6i(x) - sin OA 2 

XCKU = { t*G ) ' (83) 

We take the masses fhi at the GUT scale, and their errors Srrii, from Table II of |21) ; those 
masses refer to the MSSM with tan f3 — 10 and a GUT scale of 2 x 10 16 GeV and have 
been obtained through the renormalization-group evolution of the masses given in j22] at 
the Z°-mass scale. As for sines of the CKM angles, sin 6i±5 sin 6i, we use Table 1 in [TTj . 
We do not take into account the CKM phase in our fitting procedure; this omission will 
be justified later. 

In order to get a better understanding of our mass matrices, we perform separate 
minimizations of Xmasses an d of Xtotai- We also test the "real" versus the "complex" case. 

For the numerical mult i- dimensional minimization of the x 2_ functions we employ the 
downhill simplex method Because the problem is highly non-linear, we expect the 
existence of many local minima. 3 We start with randomly generated initial simplices. At 
the points where the numerical algorithm stops, we iterate the procedure with random 
perturbations in order to find a lower x 2 - I n this way we can be fairly certain about the 
distribution of the local minima and about the position of the global minimum. 

In the description of the fits, the concept of "pull" with respect to an observable O is 
useful. The pull of O is defined as 

pdU(0) = £^, (84) 

where the experimental value of the observable is O ± 50, while O (x) is the theoretical 
prediction of O, given as a function of the parameter set x; x is the parameter set at a 
local minimum of \ 2 ' . Thus, 

min x 2 (x) = X 2 (x) = £ [pull (O)} 2 . (85) 

o 



7 Results 

We have performed all fits and tests of our scenario separately for real and complex 
coupling matrices S and A — see Section |21 It turns out that there are no significant 
numerical differences between the two cases. The extra two phases in the "complex" case 
are unable to significantly improve our fits. Therefore, for simplicity in the following we 
confine ourselves to the "real" case. 

3 The concept "local minimum" is not understood in a strict mathematical sense, rather it refers to a 
point where the minimization algorithm successfully stops. 
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5 10 15 20 25 30 35 40 

2 

Figure 1: The distribution of local minima of Xmasses ~ Xd + x\ in the Xd~Xu plane (we 
force xi to be always negligibly small). The straight lines refer to constant values of 
Xmasses- This figure refers to the "real" case. 

Fits of the masses alone: Firstly we omit the CKM angles and test whether, with 
the mass matrices Q, we are able to fit the charged- fermion mass values at the SUSY 
GUT scale given in j2l] . In Fig. we show the distribution, in the XdrXu plane, of the 
local minima of Xmasses f° r which min Masses — 40. Though the density of points in that 
figure depends sensitively on the number of random perturbations and on the number 
of restarts of the downhill simplex procedure, the overall picture is clear. The absolute 
minimum of Xmasses * s located at x\ — ®, Xd — 23.3; the corresponding fit masses, and 
the pulls, are given in Tabled For comparison, we also show in Tabled the central mass 
values of [21]. Looking at the pulls, we see that this mass fit fails only in the mass of the 
down quark; that particular pull is responsible for almost the complete Xmasses = 23.3. A 
glance at Fig. ^also reveals that there are local minima with xt — 25 and — 1; those 
minima, the best of which is also displayed in Table ^ give rather good fits for all the 
down-type-quark masses, but fail severely in fitting the top-quark mass: the fit value is 
about one order of magnitude, or five a, smaller than the experimental value — see Tabled 
Thus, with our scenario we cannot even fit all the charged-fermion masses. However, as 
stressed in the introduction, our scenario is extreme in that it allows only for the Yukawa 
couplings of one 10 © 120 scalar representation. If we allow for small perturbations of 
the mass matrices, there are several ways out: there could be contributions from Yukawa 
couplings of one 126 several 10 and/or 120 of scalars, radiative corrections, or non- 
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0.4112 


-4.74 


1.430 


-0.321 


m s 


29.95 


29.54 


-0.090 


29.35 


-0.132 


m b 


1063.6 


1187.0 


0.873 


1188.2 


0.882 


m u 


0.7238 


0.7249 


8 x 10~ 3 


0.7321 


0.061 


m c 


210.33 


212.47 


0.113 


214.66 


0.228 


m t 


82433 


78466 


-0.269 


8778 


-4.99 



Table 1: Results of the fit for the "real" case without the CKM angles. The values of 
the masses in the second column, and the corresponding errors Sm for the calculation of 
the pulls, have been taken from [21]. The third column gives our best fit and the fourth 
column displays the corresponding pulls. The fifth and sixth columns refer to best fit in 
the region of small x% the region in Fig. ^with \\ — 1 an d x\ — 25. All the masses 
are in units of MeV. 

renormalizable terms. Consequently, the absolute minimum of Xmasses can b e considered 
acceptable, since it fails only for m^, which is small anyway. On the other hand, our 
philosophy of small perturbations forces us to discard the local minimum where the fit 
value of m t is one order of magnitude too small. 

Fits with masses and CKM angles: Figure 121 shows the distribution in the Xd~Xu 
plane, for the "real" case, of the local minima of Xtotai which have minx 2 otal < 50. We see 
that the gross feature — the lower left corner is devoid of local minima — is the same as in 
the fit without CKM angles. The previous local minimum at x\ — 25 and — 1 is now 
the absolute minimum. That absolute minimum is given in detail in Table 121 We see that 
the fit of m t is as unacceptably bad as before, but the CKM angles are well reproduced. 

Moving to the zone of x\ ~ 1 i n Fig- 121 where the top-quark mass is well reproduced, 
we find the following characteristic features: 

• In that zone the best fit has Xtotai — 45.4. The corresponding fit values and pulls 
are shown in columns five and six of Table |2l respectively. 

• Varying Xtotai between 45.4 and 49, we find that the pull of changes roughly 
from —2 to —3. 

• For that range of Xtotab ^ ne P u ^ °f m s remains close to —6, i.e. the fit value of m s 
is one order of magnitude lower than the experimental value — indeed, m s turns out 
hardly larger than m^. This is the main reason why Xtotai i s so i n the region of 
low xl- 
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Figure 2: The distribution of local minima of Xtotai ( m the "real" case) in the Xd~Xu plane. 
The straight lines refer to constant Xa + Xu — Xtotai- 

• The pull of rrib is always about +1. 

Since we are unable to reproduce well all the quark masses, we cannot expect to obtain 
a realistic CKM phase, and we have not included it in our fit. 

In view of our philosophy, we have also tried a fit of the masses and of the CKM mixing 
angles while allowing for artificially large errors in the light-fermion masses. Taking for 
instance Srrii = 5MeV for i = e, d, u, we are able to achieve Xtotai — 25.4. This is not 
really an improvement when compared to the best fit in columns 3 and 4 of Table 121 
However, the characteristics of this fit are different from those of that best fit: the pulls 
of rrib and of m t are approximately +1 and —1, respectively, whereas the fit value of m s is 
4.65 & too low. Thus, the fit with drastically increased errors in the light-fermion masses 
rather resembles the fit of columns 5 and 6 of Table El 

A numerical test of b—r unification: In Section 0] we have traced b-r unification 
to some inequalities involving the charged-lepton and the down-type-quark masses; those 
inequalities are conditions on the masses necessary for Eqs. (jUJ), ([TUj) . (JHJ), and (fTKJ) to 
have a solution. Neglecting the masses of the first generation, i.e. setting ~ m e ~ 0, 
those conditions are reformulated as [TH] 



1 - 



m u + m, 




«V + m s 



(86) 
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Xtota' 


= 26.9 


Xtotal 


= 45.4 




m 


m (x) 


pull (m) 


m (x) 


pull (m) 


m e 


0.3585 


0.3585 


-1 x 10~ 3 


0.3585 


7 x 10~ 3 


«v 


75.67 


75.67 


2 x 10~ 3 


75.67 


4 x 10" 3 


m T 


1292.2 


1292.2 


-4 x 10- 3 


1292.2 


-3 x 10" 3 


m d 


1.504 


1.563 


0.141 


1.044 


-1.993 


m s 


29.95 


28.24 


-0.376 


1.36 


-6.29 


m b 


1063.6 


1191.0 


0.903 


1225.4 


1.145 


m u 


0.7238 


0.7243 


3 x 10~ 3 


0.7279 


0.030 


m c 


210.33 


215.35 


0.264 


216.83 


0.342 


m t 


82433 


8179 


-5.03 


74145 


-0.56 




sin 9 


sin 9 (x) 


pull (sin 9) 


sin 9 (x) 


pull (sin 9) 


sin 9 12 


0.2243 


0.2242 


-0.047 


0.2243 


2 x 10" 4 


sin #23 


0.0351 


0.0348 


-0.208 


0.0352 


0.093 


sin 


0.0032 


0.0036 


0.740 


0.0034 


0.318 



Table 2: Results of the fit for the "real" case including the CKM angles. The best fit is 
described by columns three and four: this point lies at x\ — 1> Xu — 25. Columns five 
and six refer to the best fit in the region of small \\. 



Using the values for the masses in [21], this reads 

0.92 < — < 1.08. (87) 
m T 

We have performed a \ 2 analysis to check the inequalities (jHHJ). For this purpose, we 
consider 

X*(x,m h )= E 1 5mi ) + ( o,qi W6 1 • ( 88 ) 



i=e,/j,,T,d. 

This x 2 "f unc ti° n is identical with Xd + xj apart from the term corresponding to the 
bottom-quark mass, wherein we leave m b free and assign to it a very small error bar of 
1%. Then we define a minimal \ 2 as a function of m b . 



Xm (m b ) = mm X m (x, m b ) . (89) 



This function allows one to test the down-type- quark and charged-lepton mass fits with 
respect to variations of m b |21j. In Fig.0]we have plotted \di ( m b) against m b /m T ] we have 
used the mean value m T = 1292.2 MeV given in [21]. We see that exactly in the range of 
Eq. ()87|) the minimum of x\i i m b) is, for all practical purposes, zero. This confirms our 
analytic derivation of b-r unification. 

We can also test Eq. ()46)1 against our numerics. We find that that equation is repro- 
duced fairly well whenever the fit of m t is good. 
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0.80 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 

mJm 

b X 



Figure 3: The minimum of the x 2 -fu n ction of Eq. (|KH|) as a function of m,b- 

8 Summary 

In this paper we have investigated a SUSY SO (10) scenario in which the charged-fermion 
masses are generated exclusively by the renormalizable Yukawa couplings of the fermions 
to one representation 10 © 120 of scalars. We have studied the three-generations case, 
confirming the b-r unification which had previously been proved for two generations [18 . 
However, our tests of this scenario against the charged-fermion masses and against the 
CKM mixing angles at the GUT scale show that it is not satisfactory: the fit value of m t 
comes out much too low for the best fit; allowing for a larger Xtotab we are & ble to obtain 
a good fit of m t , at the price of m s turning out one order of magnitude too low and of wid 
also being too small. We thus find that the scenario investigated here is too restrictive: 
an additional mechanism for charged-fermion mass generation is required. 
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a Ciencia e a Tecnologia through the project U777-Plurianual. 
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